Search for: All records

ABSTRACT We study the morphology of hundreds of simulated central galaxies in the stellar mass range $$M_\star =$$ 107.5–1011  $$\rm M_\odot$$ from the FIREbox cosmological volume. We demonstrate that FIREbox is able to predict a wide variety of morphologies, spanning from disc-dominated objects to spheroidal galaxies supported by stellar velocity dispersion. However, the simulations predict a strong relation between morphology (degree of rotational support) and stellar mass: galaxies comparable to the Milky Way are often disc-dominated while the presence of stellar discs mostly vanishes for dwarfs with $$M_\star < 10^9 ~$$\rm M_\odot$$. This defines a ‘morphology transition’ regime for galaxies with $$10^9 < M_\star /\rm {M_\odot }< 10^{10}$$ in which discs become increasingly common, but below which discs are rare. We show that burstiness in the star formation history and the deepening of the gravitational potential strongly correlate in our simulations with this transition regime, with discs forming in objects with lower levels of burstiness in the last $$\sim 6$$ Gyr and haloes with mass $$\sim 10^{11} ~ \rm {{\rm M}_{\odot }}$$ and above. While observations support a transition towards thicker discs in the regime of dwarfs, our results are in partial disagreement with observations of at least some largely rotationally supported gas discs in dwarfs with $$M_\star < 10^9$$\rm M_\odot$$. This study highlights dwarf morphology as a fundamental benchmark for testing future galaxy formation models. 

Context.Disc-dominated galaxies can be difficult to accommodate in a hierarchical formation scenario such as Λ cold dark matter (ΛCDM), where mergers are an important growth mechanism. However, observational evidence indicates that these galaxies are common in the Universe. Aims.We seek to characterise the conditions that lead to the formation of disc-dominated galaxies within ΛCDM. Methods.We used dynamical decomposition of the stellar particles in all galaxies with stellar massM_∗= [10¹⁰− 10¹¹] M_⊙within the cosmological hydrodynamical simulation Illustris TNG100. We selected a sample of 43 mostly-disc galaxies that have less than ∼10% of their mass in a bulge component. For comparison, we also studied two additional stellar-mass matched samples: 43 intermediate galaxies having ∼30% of their stellar mass in the bulge and 43 with a purely spheroidal-like morphology. Results.We find that the selection purely based on stellar dynamics is able to reproduce the expected stellar population trends of different morphological types, with higher star-formation rates and younger stars in disc-dominated galaxies. Halo spin seems to play no role in the morphology of the galaxies, in agreement with previous works. At a fixedM_*, our mostly-disc and intermediate samples form in dark matter haloes that are two to ten times less massive than the spheroidal sample, highlighting a higher efficiency in disc galaxies to retain and condensate their baryons. On average, mergers are less prevalent in the buildup of discs than in spheroidal galaxies, but there is a large scatter, including the existence of mostly-disc galaxies, with 15%–30% of their stars coming from accreted origin. Discs start to form early on, settling their low vertical velocity dispersion as early as 9–10 Gyr ago, although the dominance of the disc over the spheroid was established more recently (3–4 Gyr lookback time). The most rotationally supported discs form in haloes with the lowest virial mass in the sample and the best aligned distribution of angular momentum in the gas. 

Abstract Field dwarf galaxies not actively forming stars are relatively rare in the local Universe, but are present in cosmological hydrodynamical simulations. We use the TNG50 simulation to investigate their origin and find that they all result from environmental effects that have removed or reduced their gas content. Quenched field dwarfs consist of either backsplash objects ejected from a massive host or of systems that have lost their gas after crossing overdense regions such as filaments or sheets (“cosmic web stripping”). Quenched fractions rise steeply with decreasing stellar mass, with quenched systems making up roughly ∼15% of all field dwarfs (i.e., excluding satellites) with stellar masses 10⁷ < M_⋆/M_⊙ < 10⁹. This fraction drops to only ∼1% when a strict isolation criterion that requires no neighbors withM_⋆ > 10⁹M_⊙within 1.5 Mpc is applied. Of these isolated dwarfs, ∼6% are backsplash, while the other ∼94% have been affected by the cosmic web. Backsplash systems are more deficient in dark matter, have retained less or no gas, and have stopped forming stars earlier than cosmic web-stripped systems. The discovery of deeply isolated dwarf galaxies that were quenched relatively recently would lend observational support to the prediction that the cosmic web is capable of inducing the cessation of star formation in dwarfs. 

Abstract We study the stellar properties of a sample of simulated ultradiffuse galaxies (UDGs) with stellar massM_⋆= 10^7.5–10⁹M_⊙, selected from the TNG50 simulation, where UDGs form mainly in high-spin dwarf-mass halos. We divide our sample into star-forming and quenched UDGs, finding good agreement with the stellar assembly history measured in observations. Star-forming UDGs and quenched UDGs withM_⋆≥ 10⁸M_⊙in our sample are particularly inefficient at forming stars, having 2–10 times less stellar mass than non-UDGs for the same virial mass halo. These results are consistent with recent mass inferences in UDG samples and suggest that the most inefficient UDGs arise from a late assembly of the dark matter mass followed by a stellar growth that is comparatively slower (for star-forming UDGs) or that was interrupted due to environmental removal of the gas (for quenched UDGs). Regardless of efficiency, UDGs are 60% poorer in [Fe/H] than the population of non-UDGs at a fixed stellar mass, with the most extreme objects having metal content consistent with the simulated mass–metallicity relation atz∼ 2. Quenched UDGs stop their star formation in shorter timescales than non-UDGs of similar mass and are, as a consequence, alpha enhanced with respect to non-UDGs. We identify metallicity profiles in UDGs as a potential avenue to distinguish between different formation paths for these galaxies, where gentle formation as a result of high-spin halos would present well-defined declining metallicity radial profiles while powerful-outflows or tidal stripping formation models would lead to flatter or constant metallicity as a function of radius due to the inherent mixing of stellar orbits. 

Abstract Simulations of galaxy formation are mostly unable to resolve the energy-conserving phase of individual supernova events, having to resort to subgrid models to distribute the energy and momentum resulting from stellar feedback. However, the properties of these simulated galaxies, including the morphology, stellar mass formed, and the burstiness of the star formation history, are highly sensitive to the numerical choices adopted in these subgrid models. Using the SMUGGLE stellar feedback model, we carry out idealized simulations of anM_vir∼ 10¹⁰M_⊙dwarf galaxy, a regime where most simulation codes predict significant burstiness in star formation, resulting in strong gas flows that lead to the formation of dark matter cores. We find that by varying only the directional distribution of momentum imparted from supernovae to the surrounding gas, while holding the total momentum per supernova constant, bursty star formation may be amplified or completely suppressed, and the total stellar mass formed can vary by as much as a factor of ∼3. In particular, when momentum is primarily directed perpendicular to the gas disk, less bursty and lower overall star formation rates result, yielding less gas turbulence, more disky morphologies, and a retention of cuspy dark matter density profiles. An improved understanding of the nonlinear coupling of stellar feedback into inhomogeneous gaseous media is thus needed to make robust predictions for stellar morphologies and dark matter core formation in dwarfs independent of uncertain numerical choices in the baryonic treatment. 

ABSTRACT We use the TNG50 from the IllustrisTNG suite of cosmological hydrodynamical simulation, complemented by a catalogue of tagged globular clusters, to investigate the properties and build up of two extended luminous components: the intra-cluster light (ICL) and the intra-cluster globular clusters (ICGCs). We select the 39 most massive groups and clusters in the box, spanning the range of virial masses $$5 \times 10^{12} \lt \rm M_{200}/\rm {\rm M}_{\odot } \lt 2 \times 10^{14}$$. We find good agreement between predictions from the simulations and current observational estimates of the fraction of mass in the ICL and its radial extension. The stellar mass of the ICL is only $$\sim 10~{{\ \rm per\ cent}}$$–20 per cent of the stellar mass in the central galaxy but encodes useful information on the assembly history of the group or cluster. About half the ICL in all our systems is brought in by galaxies in a narrow stellar mass range, M* = 1010–1011 M⊙. However, the contribution of low-mass galaxies (M* < 1010 M⊙) to the build up of the ICL varies broadly from system to system, $$\sim 5~{{\ \rm per\ cent}}-45~{{\ \rm per\ cent}}$$, a feature that might be recovered from the observable properties of the ICL at z = 0. At fixed virial mass, systems where the accretion of dwarf galaxies plays an important role have shallower metallicity profiles, less metal content, and a lower stellar mass in the ICL than systems where the main contributors are more massive galaxies. We show that intra-cluster GCs are also good tracers of this history, representing a valuable alternative when diffuse light is not detectable. 

We study the fraction of the intra-cluster light (ICL) formed in-situ in the three most massive clusters of the TNG50 simulation, with virial masses $\sim {10}^{14}$. We find that a significant fraction of ICL stars ( $8\%$- $28\%$) are born in-situ. This amounts to a total stellar mass comparable to the central galaxy itself. Contrary to simple expectations, only a sub-dominant fraction of these in-situ ICL stars are born in the central regions and later re-distributed to more energetic orbits during mergers. Instead, many in-situ ICL stars form directly hundreds of kiloparsecs away from the central galaxy, in clouds condensing out of the circum-cluster medium. The simulations predict a present-date diffuse star formation rate of $$1 ${\mathrm{M}}_{\odot }$/yr, with higher rates at higher redshifts. The diffuse star forming component of the ICL is filamentary in nature, extends for hundreds of kiloparsecs and traces the distribution of neutral gas in the cluster host halo. We discuss briefly how numerical details of the baryonic treatment in the simulation, in particular the density threshold for star formation and the equation of state, may play a role in this result. We conclude that a sensitivity of $1.6\times {10}^{-19}-2.6\times {10}^{-18}$erg s ${}^{-1}$cm ${}^{-2}$arcsec ${}^{-2}$in H ${}_{\alpha }$flux (beyond current observational capabilities) would be necessary to detect this diffuse star-forming component in galaxy clusters. 

ABSTRACT Stellar feedback plays a crucial role in regulating baryon cycles of a galactic ecosystem, and may manifest itself in the formation of superbubbles in the interstellar medium. In this work, we used a set of high-resolution simulations to systematically study the properties and evolution of superbubbles in galactic environments. The simulations were based on the SMUGGLE galaxy formation framework using the hydrodynamical moving-mesh code arepo, reaching a spatial resolution of $$\sim 4 \, \rm pc$$ and mass resolution of $$\sim 10^3 \, \rm M_{\odot }$$. We identified superbubbles and tracked their time evolution using the parent stellar associations within the bubbles. The X-ray luminosity-size distribution of superbubbles in the fiducial run is largely consistent with the observations of nearby galaxies. The size of superbubbles shows a double-peaked distribution, with the peaks attributed to early feedback (radiative and stellar wind feedback) and supernova feedback. The early feedback tends to suppress the subsequent supernova feedback, and it is strongly influenced by star formation efficiency, which regulates the environmental density. Our results show that the volume filling factor of hot gas (T > 105.5 K) is about $$12~{{\ \rm per\ cent}}$$ averaged over a region of 4 kpc in height and 20 kpc in radius centred on the disc of the galaxy. Overall, the properties of superbubbles are sensitive to the choice of subgrid galaxy formation models and can, therefore, be used to constrain these models. 

ABSTRACT In this study, we modify the semi-analytic model galacticus in order to accurately reproduce the observed properties of dwarf galaxies in the Milky Way. We find that reproducing observational determinations of the halo occupation fraction and mass–metallicity relation for dwarf galaxies requires us to include H2 cooling, an updated ultraviolet background radiation model, and to introduce a model for the metal content of the intergalactic medium. By fine-tuning various model parameters and incorporating empirical constraints, we have tailored the model to match the statistical properties of Milky Way dwarf galaxies, such as their luminosity function and size–mass relation. We have validated our modified semi-analytic framework by undertaking a comparative analysis of the resulting galaxy–halo connection. We predict a total of $$300 ^{+75} _{-99}$$ satellites with an absolute V-band magnitude (MV) less than 0 within 300 kpc from our Milky Way analogues. The fraction of subhaloes that host a galaxy at least this bright drops to 50 per cent by a halo peak mass of ∼8.9 × 107 M⊙, consistent with the occupation fraction inferred from the latest observations of Milky Way satellite population. 

ABSTRACT The velocity dispersion of globular clusters (GCs) around ultra-diffuse galaxies (UDGs) in the Virgo cluster spans a wide range, including cases where GC kinematics suggest haloes as massive as (or even more massive than) that of the Milky Way around these faint dwarfs. We analyse the catalogues of GCs derived in post-processing from the TNG50 cosmological simulation to study the GC system kinematics and abundance of simulated UDGs in galaxy groups and clusters. UDGs in this simulation reside exclusively in dwarf-mass haloes with M200 ≲ 1011.2 M⊙. When considering only GCs gravitationally bound to simulated UDGs, we find GCs properties that overlap well with several observational measurements for UDGs. In particular, no bias towards overly massive haloes is inferred from the study of bound GCs, confirming that GCs are good tracers of UDG halo mass. However, we find that contamination by intracluster GCs may, in some cases, substantially increase velocity dispersion estimates when performing projected mock observations of our sample. We caution that targets with less than 10 GC tracers are particularly prone to severe uncertainties. Measuring the stellar kinematics of the host galaxy should help confirm the unusually massive haloes suggested by GC kinematics around some UDGs.